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^ . Abstract 

00 
00 

^^ ■ 'current' u, d quarks interacting pairwise via gluon-exchange- propagator in its non- 

^^ , perturbative regime, is employed for the calculation of baryon self-energy via quark- 

^^ I loop integrals. To that end the baryon-ggg vertex function is derived under Covari- 

ant Instantaneity Ansatz (CIA), using Green's function techniques. This is a 3-body 
extension of an earlier qq (2-body) result on the exact 3D-4D interconnection for 
the respective BS wave functions under 3D kernel support, precalibrated to both qq 
^j' and qqq spectra plus other observables. The quark loop integrals for the neutron (n) 

r^ • - proton (p) mass difference receive contributions from : i) the strong SU(2) effect 

arising from the d — u mass difference (4 MeV); ii) the e.m. effect of the respective 
quark charges. The resultant n — p difference comes dominantly from d — u effect 
{+1.71Mev), which is mildly offset by e.m. effect (—0.44), subject to gauge correc- 
tions. To that end, a general method for QED gauge corrections to an arbitrary 
momentum dependent vertex function is outlined, and on on a proportionate basis 
from the (two-body) kaon case, the net n-p difference works out at just above 1 
MeV. A critical comparison is given with QCD sum rules results. 
PACS : 11.10 st ; 12.35 Ht ; 12.70 -|- q Key-words: baryon self-energy; qqq BS-vertex 
fn; Cov. Inst. Ansatz; 3D-4D interconnection; Green's fn method; QED gauge inv; 
n-p diff. 
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1 Introduction: Relativist ic 2- and 3-Quark Hadrons 

Soon after the advent of the Faddeev theory [1], the relativistic 3-body problem [2] at- 
tracted instant attention as a non-trivial dynamical problem, as distinct from earlier 
"kinematical" attempts [3] at a relativistic formulation of its wave function. In this re- 
spect the relativistic 3-baryon problem had been more of academic than practical interest 
(until the 'pion' got involved as a key ingredient), but the situation changed qualitatively 
when this 3-body problem started being viewed at the quark level. Looking back after 
25 years it appears that the first serious attempt in this direction was made by Feynman 
et al [4] who gave a unified formulation of both the qq (meson) and qqq (baryon) prob- 
lems under a common dynamical framework, bringing out rather sharply an underlying 
duality between these two systems which in turn signifies a more basic duality between 
a qq diquark [5] and a q antiquark. Indeed the diquark description is quite compact and 
adequate for many practical purposes involving the baryon, but the more microscopic 
qqq description which brings out the fuller permutation {S3) symmetry in the baryon is 
necessary for the actual details of a full-fledged dynamical treatment [4]. 

1.1 BS Dynamics: 3D vs 4D Forms 

Although the FKR theory [4] marked the first step in this direction, it suffered from an 
inadequate treatment of the time-like d.o.f. which showed up in several ways. The lat- 
ter has by itself a long history of attempts at 3D formulations of the 4D Bethe-Salpeter 
Equation (BSE) for qq qqq systems:Instantaneous Approximation (lA) [6a]; Quasipoten- 
tials [6b]; On-shellness of the associated propagators [6c]; and others [6d]. Some of these 
approaches have been reviewed elsewhere [7] in the context of a unified BS formulation of 
qq and qqq in the FKR [4] spirit. 

A different form of 3D reduction, which is of more recent origin [8,9] is based on 
the Markov- Yukawa Transversality Condition [10] in which the BSE kernel is given a 3D 
support by demanding that it be a function of the relative momentum q transverse to the 
total 4-momentum P. This condition, termed the "Covariant Instaneity Ansatz" (CIA) 
[8], is somewhat complementary to the more conventional approaches [6] in which the 
BS kernel is left untouched but the propagators are manipulated in various ways to give 
a 3D reduction to the BSE, whereas in the Markov- Yukawa [8-10] method, which has 
been termed CIA [8], it is the other way around. In the conventional 3D reductions [6] 
there is no going back to the original 4D form, whereas in the CIA [8], the two forms are 
fully interchangeable according to need, thus offering a possible Lorentz covariant way to 
reconcile the apparently conflicting demands of spectroscopy [11] needing a 3D BSE, with 
the 4D BS vertex functions needed for quark-loop integrals. The effectiveness of the CIA 
in giving a concrete shape to such a "two-tier" philosophy of spectra-cum-loop integrals 
was summarised in a semi-review [12] in the form of appropriate BSE's for both qq and 
qqq systems with vector- type kernels [7] with 3D support, albeit with slight modifications 
[7] in the respective BSE structures to facilitate greater 'manoeuvreability', in the spirit 
of similar efforts [13] in the past. Further, the observed spectroscopy [11] is well satisfied 
on both qq [14] and qqq [15] sectors with a common set of parameters for the respective 
kernels (the qq kernel has just half the strength of the qq kernel due to color effects), so 
that the respective vertex functions are entirely determined within the CIA formalism. 



1.2 QCD-Motivated BSE in 3D-4D Form 

The other aspect of this 'two-tier' formahsm concerns the crucial property of chiral sym- 
metry and its dynamical breaking. The first part (chiral symmetry) is ensured without 
extra charge by the vector character of the kernel that had been present all along in this 
program [7,16], since the BS-kernel is a direct reflection of an effective 4-fermion term in 
the input Lagrangian. Indeed the vector type character of the latter lends a natural gluon 
exchange flavour to such a pairwise interaction among 'current' (almost massless) u, d 
quarks at the Lagrangian level. This structure is quite general [17], and can be adapted 
to the QCD requirements on the gluonic propagator involved in the pairwise interaction 
kernel. Of this, the perturbative part (which is well understood) is quite explicit, but the 
non-perturbative (infrared) part is not yet derivable from formal QCD premises [18]. It 
can nevertheless be simulated in a sufficiently realistic manner at the phenomenological 
level [7,19], so as to satisfy the standard constraints of confinement as well as explicit 
QCD features [20] in terms of a basically 3D BSE kernel structure. 

The second part, viz., dynamical breaking of chiral symmetry (DBxS) is implemented 
via the Nambu-Jonalasino mechanism [21] whose full-fiedged form amounts to adopting 
the 'non-trivial' solution of the Schwinger- Dyson Equation (SDE) derived from a given 
input, chirally symmetric Lagrangian with current quarks. A mass function m{p) [17,18] 
is thus generated whose low-momentum value may be identified with the bulk of the 
'constituent' mass (m^) of the u, d quarks. This accords with Politzer additivity [22], viz., 
rriq = m{0)+mc; where rric, the current mass, is small. This was also shown in the context 
of a BSE-cum-SDE treatment [23] within the CIA formalism [8]. Thus formally the BS- 
kernel may be regarded as a non-perturbative gluon propagator [23] in a BSE framework 
involving the dynamical/constituent mass [19, 21-24] in the quark propagator. 

To recapitulate, the CIA which gives an exact interconnection between the 3D and 
4D forms of the BSE, provides a unified view of 2- and 3-quark hadrons, its 3D reduc- 
tion being meant for spectroscopy [14-15], and the reconstructed 4D form [12,13] for 
identifying the respective hadron-quark vertex functions as the key ingredients for 4D 
quark-loop integrals. The formalism stems from a strongly QCD-motivated Lagrangian 
with current quarks whose pairwise interaction is mediated by a gluonic propagator in 
its non-perturbative regime. The QCD feature of chiral symmetry is ensured by the 
vector nature of this interaction, while its dynamical breaking is the result of a non- 
trivial solution of the SDE [17,23]. Thus, unlike in conventional potential models [25], 
the constituent mass so generated is not a phenomenological artefact, but the result of a 
self-consistent solution of the SDE [17, 19, 23], so that the standard (constituent) mass 
employed for spectroscopy [14-15] 'checks' with the output dynamical mass at low mo- 
mentum [23]. Thus there are only two genuine input parameters Co, ujq, that characterize 
the (phenomenological) structure of the non-perturbative gluon propagator which serves 
for both the 2- and 3-quark spectra in a unified fashion [14-15]. In this formalism, these 
two constants play a role somewhat similar to that of the (input) 'condensates' in the 
theory of QCD sum rules [26]. 

1.3 Comparison With Chiral Perturbation Theory, Etc 

Before proceeding further, let us pause to compare this approach with other dynamical 
methods, e.g., chiral perturbation theory [27] which has more explicit QCD features, albeit 



in the perturbative regime, leading to expansions in the momenta. This is a powerful 
theoretical approach employing the (chiral) symmetry of QCD; its essential parameters 
are the current quark masses, and the method works very efficiently where its premises 
are logically applicable. Thus it predicts the ground state spectra of light quark hadrons, 
including their mass splittings due to strong and e.m. breaking of SU(2), but not the 
spectra of L-excited hadrons. The latter on the other hand demand a "closed form" 
approach to incorporate the "soft" off-shell effects which in turn require a non-trivial 
handle on the infrared (non-perturbative) part of the gluonic propagator, something which 
the present state of the QCD art does not yet provide. Thus one needs a phenomenological 
input even in standard BSE-SDE approaches [18], as discussed elsewhere [23]. The chiral 
perturbation theory [27] also lacks this vital ingredient, as seen from the absence of form 
factors in its 'point' Lagrangians [27] with at most derivative terms. This shows up, e.g., 
through its inability to predict L-excited spectra, and finer aspects (such as convergence) 
of 4D quark-loop integrals which depend crucially on these "off-shell" features. Physically 
this amounts to the absence of a 'confinement scale' which governs these form factors. In 
other BSE-cum-SDE approaches [17-19], including the present 'two-tier' CIA formalism 
[8,12], this 'scale' is an integral part of the structure of the non-perturbative part of 
the gluon propagator [19,23], with a built-in QCD feature of chiral symmetry and its 
dynamical breaking through the non-trivial solution of the SDE [17,19,23]. This not only 
facilitates the prediction of L-excited spectra [19,14-15] but also provides a form factor 
for the hadron-quark vertex function which greatly enhances its applicability to various 
4D quark-loop integrals; see [8, 23-24, 28]. 

1.4 Application to n-p Mass Difference with 3D-4D qqq BSE 

After this excursion on the philosophy of this two-tier BSE approach, vis-a-vis some others 
[26,27], we may now state the objective of the present paper: A typical application of the 
4D baryon-ggg wave function reconstructed [12] from the 3D qqq BSE, to the n — p mass 
difference, as a 3-body generalization of the corresponding qq-meson problem [28]. Unlike 
the 2-body case, however, where the 3D-4D interconnection is exactly reversible [8], a 
4D reconstruction for a 3-body system involves a loss of information on the 4D Hilbert 
space, so that the reversal of steps is in principle non unique, and requires a ID 5-function 
to fill up the information gap between 3D and 4D Hilbert space which may be directly 
attributed to the CIA ansatz of a 3D support to the pairwise kernel. The 2-body case 
just escapes this pathology as it represents a sort of degenerate situation, but the price of 
a 3D kernel support must show up in a reconstruction of the 4D BSE from its reduced 3D 
form in any (n > 2)-body problem [29]. A plausible 'CIA' structure for the 4D qqq wave 
function was suggested in [12] in a semi- intuitive fashion, but a more formal mathematical 
basis has since been found [29] through the use of Green's function techniques, so that 
the reconstructed 4D form reduces exactly to the (known) 3D form as a consistency check 
[29]. The final result, which is almost the same as the earlier conjecture, eq.(5.15) of [12], 
except for a constant that does not affect the normalization, contains a ID 5-function 
corresponding to the on-shell propagation of the spectator between two successive vertex 
points. As explained in detail in [29], this ID 6- function must not be confused with any 
signature of "non-connectedness" in the 3-body wave function [30], since the 3D form 
is fully connected. Rather, it can be likened to a (Fermi-type) (5-function potential in 
estimating the effect of chemical binding on the scattering of very slow neutrons by a 



hydrogen molecule [31], as a practical devise to fill up the vast mismatch of scales in the 
interactions at nuclear vs molecular levels. In any case the ID (5-function appearing in 
this structure is entirely innocuous as it gets integrated out in any physical (quark loop) 
amplitude including the BS normalization (see Sec. 2 below). 

Now to recall the physics of the n-p mass difference, this quantity receives contributions 
of opposite signs from two main sources: i) a positive one from the strong SU(2) d — u 
mass difference; ii) a negative one from e.m. splittings. A third source, the effect of quark 
condensates, which plays a crucial role in QCD sum rule studies [26] gives rather small 
contributions in this non-perturbative approach, as already found in the meson case [28], 
and will therefore be neglected; for a detailed discussion on this issue, see sec. 5. 2. 

On the other hand the problem of gauge invariance (g.i.) of the e.m. contribution (ii) 
is a more tricky issue in view of the otherwise arbitrary nature of the extended form factor 
associated with the baryon — qqq vertex function, in contrast with the point-like vertices 
involved, e.g., in the corresponding QCD-sum rule studies [32]. For the same reason, the 
g.i. issue could not be addressed in [28] for the meson case, in the hope that it would 
be of the same order of magnitude as the relatively small (20of the e.m. contribution 
[28]. The g.i. problem with arbitrary hadron-quark vertex functions is best left to a 
separate, more substantial investigation. In the meantime in this paper we shall estimate 
the g.i. corrections for a 2-body problem on the lines of [32] adapted to an arbitrary 
meson- quark vertex function. The steps are indicated in an appendix (Appendix C) for 
the kaon problem as a test case, the results of which are provisionally considered as an 
indication of the nature of the g.i. corrections to be expected for the 3-body problem on 
hand, pending a formal treatment later. The kaon result indicates an increase of 0.612 
MeV in the e.m. self-energy (1.032 MeV) arising from fig 1(b) of [28]. 

1.5 Contents of the Paper 

In Sec. 2 we collect the various pieces of the central quantity of the present investigation, 
viz., the 4D baryon— qqq vertex function in terms of 3D quantities, with the inclusion of the 
spin and isospin d.o.f. on the lines of an earlier study [33]. Thus equipped, we outline the 
main steps leading to an explicit evaluation of the normalization integral, using Feynman 
diagrams shown in figs.l(a.b,c). A complex basis [9, 34, 35] for 3D momentum variables 
facilitates the evaluation of the resulting 3D x 3D integrals, after the time-like momenta 
have been eliminated by 'pole' integrations on identical lines to the corresponding qq 
problem [12,23,24]. In Sec. 3 we evaluate the 'shift' in the nucleon mass due to strong 
SU(2) breaking, by inserting a mass shift operator —SmT^^'^^/2 in place of ^7/^64 at each 
of the corresponding 7- vertices of figs.l(a,b,c), as shown in figs. 2.(a,b,c). Here Sm = 4 
MeV is the 'standard' d-u mass difference [24,28] taken as the basic input. Sec. 4 sketches 
the evaluation of the e.m. contribution in accordance with the diagrams of fig.3(a,b,c). 
The details of the e.m. approximations employed are collected in Appendix A. Appendix 
B sketches the main steps of the derivation [29] for the 4D structure, eq.(5.15) of [12], 
of the baryon-ggg vertex function by the Green's function method for 3 spinless quarks. 
Finally Appendix C sketckes the steps leading to the g.i. corrections [32] for the kaon 
case, as a sort of facsimile of similar corrections expected for the n — p problem on hand. 
Sec. 5 summarises our findings and conclusions vis-a-vis other methods. 



2 Normalization of the Baryon-qqq Vertex Function 

To outline the structure of the baryon-ggg vertex function from a CIA-governed BSE 
[12-13], we shall generally follow the notation, normalization and phase convention for 
the various symbols as given in [13], but adapted to the equal mass kinematics {mi = 
1712 = iTT-s = iTT-q)- The SU(2) mass difference 6m (^ AMeV) between d and u quarks will 
be taken into account only through a 2-point vertex [—Sm.T3^''^/2] inserted in the quark 
propagators in figs. 2 (in place of i'y^Ci for a photon), but not in the structure of the vertex 
function. The vertex function is written in three pieces in each of which one quark plays 
the role of the 'spectator' by turn. For the spin structure (not given in [13]) we employ 
the convention of [3] which was extended in [33] to incorporate the S'3-symmetry for the 
spin-cum- isospin structure in the Verde [36] notation [37]. The full 4D BS wave function 
^ reads as [13,33,34] : 



^AiAsAa = (Fi + Fa + F3) x [xV + x"0"]/V2; (2.1) 

Ai = m/ + p,2; (^ = 1,2,3). (2.2) 

Here x' ^-^id x" ^-^e the relativistic "spin" wave functions in a 2-component mixed sym- 
metric 5*3 basis which for a 56 baryon go with the associated isospin functions 0' and 0" 
respectively. These are given by [3,33] : 

[X'W,^ = [{M - i^.P)i^,C/V2]p, X U{P)J{2M) (2.3) 

[x"]/3,;a = [(M - t^.P)^,C/VQ]p, X t^,^^U{P)J{2M) (2.4) 

in a spinorial basis [3,33] in which the index a refers to the 'active' quark (interacting 
with an external photon line, fig.l), while /5, 7 characterize the other two, with the further 
convention that 7 refers to the "spectator" in a given diagram, fig.(l). The 'hat' on 
7 signifies its perpendicularity to P^, viz., 7.P = 0. The notations in eqs. (2.3-4) are 
standard, with a common Dirac basis for the entire structure, and 'C is the charge 
conjugation operator for quark 7^3 in a 23-grouping [3,33]. P^ is the baryon 4-momentum, 
U{P) is its spinor representation, and (M — i'j.P)/{2M) its energy projection operator 
[3,33]. Further, because of the full S'3-symmetry of the last factor in (2.1), the (1,2,3) 
indices can be permuted as needed for the diagram on hand. Thus in fig. la, #1(q;) 
interacts with the photon ; #2(/3) is the quark which has had a 'last' gg-interaction with 
#1(0;) before emerging from the hadronic 'blob', while #3(7) is the spectator [33]. In 
fig. lb, the roles of #1 and #2 are reversed so that, of the two 'active quarks' #1 and 
#2, #2(0;) now interacts with the photon, #l(/9) has had the last qq- interaction with 
#2(0;), while #3(7) still remains the 'spectator'. These roles are cyclically permuted, 
with two more such pairs of diagrams, fig.lc), to give an identical chance to each of the 
quarks in turn [33]. Thus there are 3 such pairs of diagrams, of which only one pair is 
shown. An identical consideration applies to figs.2(a,b) with i'-f^Ci replaced by {—5mT^/2) 
consistently. The spatial vertex functions Fj are given for i = 3 by [12] : 



F3 = NB[Di2<P/2ni] X ^[27r5(A3).A3] (2.5) 

where is the full, connected qqq wave function in 3D form, and D12 is the 3D denominator 
function of the (12) subsystem . The second factor represents the effect of the spectator 



[12] whose inverse propagator Df~^{pz) off the mass shell is just A3, eq.(2.2). The 
main steps leading to this unorthodox structure which has been derived recently via 
the techniques of Green's functions [29], are sketched for completeness in Appendix B. 
As already noted in Seel, and again explained in Appendix B, its peculiar singularity 
structure in the form of a "square-root" of a ID (5-function stems from the CIA ansatz 
of a 3D support to the pairwise interaction kernel, but it is quite harmless as the former 
will appear in a linear form in the transition amplitude corresponding to any Feynman 
diagram as in figs. 1-2. The complete expressions for D12 and are given for the equal 
mass case (with #3 as spectator) by (see [12,15]): 

D12 = Ai2(M - 0^3); A12 = 2CU12' -M\l- u^f/2 (2.6) 

^12^ = mq^ + g?2; 2gf2 =Pi -P2 (2-7) 

= ^Hpi+Pl+Pl)/2f3' = e-p/3/3^ (2.8) 

(see further below for the definition of p). 

Pt = P^^ + P^■PPJM^■, pt+pf^+P^ = (2.9) 

Ui"^ = niq^ + p. ; z/3 = Us/M^onshell) (2.10) 

The /9-parameter is defined sequentially by [14,15]: 

4 
f3* = -MuJo^as{l-mg/Mf{M- <u>); < u >^ = nig^ + 3(3^8 (2.11) 

9 

1 1 o,^^ il-mg/M)\ 

a, =as -2MCq—- — - — — , (2.12) 

M— < uj > 

67r , (M- <u>) , , 

— = 291n ^ '-; (2.13) 

AQCD = 200MeV; uo = 15SMeV; Cq = 0.29 (2.14) 

The normalization Nb, eq.(2.5), is given in accordance with the Feynman diagrams l(a,b) 
by the 4D integral ( c.f.[33]) : 



zP,/M = Y: ld%2d%-^[< 0'|(23)'(1);|0' > +- < 0"|(23)".a(1)%A;,|0" >] 

123 -^ ^^2^3 <J 

+ (1^2) (2.15) 



where the matrix element for fig. la is organized as a product of two spin-factors : a 
'23-element' expressed as a Dirac trace over the indices /9, 7; and a '1-element' (with 
suppressed index a). The associated isospin functions are shown according to (2.1). 
The contribution of fig. lb is shown symbolically by 1 <^ 2, while X]i23 indicates the sum 
over all the 3 pairs cyclically. In representing eq.(2.12) we have dropped 'cross-terms' 
like Ti*Tj^ where i 7^ j, since the presence of a v5-function in each Fj ensures that a 
simultaneous 'on-shell' energy conservation oi i ^ j spectators is not possible [33]. The 
various pieces of the matrix elements in (2.14) which can be read off from fig. la in terms 
of the spin functions (2.3-4) are as follows: 

(1):,^^ = UiP)SF{p^)^^,e^SF{p,)UiP) (2.16) 



iSp ^(p) = mg + i'j.p (2.17) 

(1)%A;m = UiP)i%l5SFiPi)il^.e,SFip,)i^,^xUiPy, (2.18) 

(23)' = Tr[C-'^5iM - i^.P)im, - i^.p2){M - i^.P)^^{m, + i^.p^)C]/%M^ (2.19) 

(23)%A = Tr[C-'%{M - ii.P){m, - ii.p2){M - i-i.P)%{m, + i-i.p^)C]/^M^ (2.20) 

The 'strength' e, of the (zero- momentum) 'photon' couphng to the quark hne Pi can be 
chosen in several ways [38]. We take here the simplest possibility, viz., Cj = 1/3 each. 
The isospin matrix element is first eliminated according to [39]: 

< (f)'\l\(f)' >=< 0" |1|0" >= 1 (2.21) 

< 0'|r3^'^l0' >= -3 < 0"r3«|0" >=< rg >(p,„) (2.22) 

Eq.(2.20) suffices for (2.4), while (2.21) will be needed for the u-d mass difference operator 
—6mT3^^^/2 ; see Sec. 3. Next, the evaluation of the traces in (2.15-18) is straightforward, 
after noting that (2.15-16), a/ter spin-averaging, are expressible as traces. The results are 

(23)'^,A = (23)%A = (m, + Mz/2)(m, + Miys)9,x (2.23) 

(1);^,, = (1)",,;^ = [2Mz/i(m, + Mz/i) + Ai]^,,P^/(MAi2) (2.24) 

where ^ is a covariant Kronecker delta w.r.t. P^, viz.. 



^uX = ^uX 



Sux - PuPx/P'; {P' = -M2) (2.25) 



Collecting all these results and simplifying we get 

Nb-' = j:Jd% ^'^\l'"'^ X J d%,D,,'<l>'[e,h + e,h] (2.26) 



123 



2i7ili = f Md(Tu[2Mui{mg + Mui) + Ai]/(MAi2A2) (2.27) 

where we have "cashed" the (5(A3)-function arising from \T^\^ against the time-like com- 
ponent of c?*^P3, and used the results 

c?^gi2 = d^uMdai2; z/i,2 = (1 - 1^3) ± ^12 (2.28) 

The integration over d(7i2 involves single and double poles arising from the propagators 
Ai 2~^ in (2.26), while the value of u^ is taken 'on-shell' at u^/M after the (5(A3)- function 
has been cashed. The result of a basic (Ti2-integration is 

f Md(Ti2Ai-^A2'^ = 2i7i/Di2 (2.29) 

from which others can be deduced by differentiation under unequal mass kinematics, or 
directly through a 'double pole' integration. The net result for Ji + 12, eq.(2.26), is given 
in eq.(2.41) below. Further, the individual terms of the summation I]i23 i^i (2.25) are 
fixed by the values chosen for Cj (which need not be specified in advance, as they can be 
adapted to other conventions too [38]; see Sec. 3). 



^z_ = i?2 e'^^ 



V2z3 = i?3 e 



ie-i 



The integration in (2.25) can be considerably simplified in a complex basis [15,34] 
defined (in momentum space) by : 

V2zi = ^i + ii]i; V2zi* = ^i- ii]i; (2.30) 

Vs^i = Pii - P2U ^Vi = -'^Psi + Pii + P2i] (2.31) 

where we now employ the alternative notation pu for p^, in view of its basically 3D 
content. In terms of Zi and Zi* , the 6D integration in (2.25) is expressed as 

d^psd^u = {V3/2fd^^d^r] = dhdh* (2.32) 

The further representation [15,34] 

d^zd'^z* = {dzj^dzJ'){dz^dz^*){dz^dz^*) (2.33) 

where 

V2z+ = Ri e'^'; V2z_* = Ri e-'^' (2.34) 

V2z+* = i?2 e-'^^ (2.35) 

V2z3* = i?3 e-'^' (2.36) 

reduces the 6D integration (2.32) merely to n^ dRi' dR2' dR^ ^ since the 6i -variables {not 
Euler angles!) are not involved in the integrands encountered, and just sum up to (27r)^. 
The positive variables i?j, (i = 1,2,3), are related to the ^i.rji variables by 

p = Rl + Rl + Rl = ^'' + rf = 2z^zC (2.37) 

To convert the variables cUj that appear in the integrals (2.28) in terms of the -Ri,2,3 
variables is a straightforward but tedious process which can be somewhat simplified in 
terms of the intermediate variables ^ — rf and 2^.77 which form a [2,1] representation [36] 
of 5'3-symmetry at the 'quadratic' level. Now because of the full 6*3 -symmetry of the 6D 
integral (2.32), together with the (fortunate) circumstance of equal mass quarks in the 
problem on hand, the integrand as a whole is S^ -symmetric which permits the following 
simplification: Each of the quantities pf and qj inside (2.32) can be expanded as 

PI2 = P/2 + ie - ^')/4 ± v^e/^/2; Pi = p/2 - (^ - v')/^ (2.38) 

ql, = 3e/^ = p/2 + ie-v')/4^ (2.39) 

In all these terms the principal quantity is p/2, while the resultant effects of the mixed- 
symmetric corrections will show up only in the fourth order, etc. In the present case of 
equal mass kinematics it is a good approximation to neglect the latter terms, as has also 
been found for the qqq mass spectral results [15], so that all quantities are expressed in 
terms of p only : 

wi,2,3 ~ UJ12 ~ ujp] Up^ = ml + p/2 (2.40) 

D^2 ^ 2{M - uj,;)[ujI - (M - uOpf/A]. (2.41) 

The rest of the integration is expressed entirely in terms of the p-variable, with the 
resultant 6D measure given by 

/ d%d%2F{p) = (7rv^/2)3 J p^dp/2F{p) (2.42) 



These considerations suffice for evaluating the integrals Ji and I2 whose resultant value 
is now given for Cj = 1/3 by : 

DUh+h) = \ml+m,{M-uj,) + {M-uj,f/A\x{M-u^f/ujp+D^2{2m,+M-ujp) (2.43) 

Substitution in (2.25) yields Nb directly. The numerical values are given coUectvely at 
the end of Sec. 4. 

3 Strong SU(2) Mass Difference for the Nucleon 

This calculation is on almost identical lines to Sec. 2, except for the substitution iei7^ 
to —SmT3^^y2 in figs.l(a,b) to give figs.2(a,b) which represent the effect of insertion of 
a 2-point vertex in a quark line. Indeed we can directly start from the counterpart of 
eq.(2.15) which gives the 'strong' mass shift as: 

^5M., = Y: fd%2d%^f^ X [< 0'|(23)'(1)'|0' > +1 < 0"|(23)%,(1)\J0" >] 

123 -^ ^^2^3 -J 

+ (1^2) (3.1) 

where we have now employed eq.(2.21) for the isospin factors, and the counterparts of 
(2.16) and (2.18) are respectively 

(1)' = f/(P)5^(pi)[-5mr3«/2]^F(pi)f/(P); (3.2) 

(1)%A = lf{P)t%-f,SF{pi)[-Smrs^'^/2]SF{p,)t-f,^xU{P) (3.3) 

while the definitions (2.19) and (2.20) remain unaltered. As a result, eq.(2.24) remains 
valid, while the counterpart of (2.23) becomes 

{lye^x = -3(1)".A = [2m,(m, + Mu,) - A,%x{-Sm/2)/Al (3.4) 

Carrying out the (iai2-integration, the result for 6Mst is now given by the counterpart of 
(2.26), viz., 

SM,, = 3iV||rf3p3^^^^|^ X j d'q,2D,2^'[J, + J2](-5mr3/6) (3.5) 

in the form of an isospin operator "ts" for the nucleon, where we have represented the 
effect of I] 123 by a factor of "3" , and 

7T1 1 

Dl^[J^ + J2] = -^[{mg + -{M-iJp)){M-u;p)[2ijl + mg{M-u;p)] 

LOp Z 

+ (M - LUpyimg + (M - iUp)/2f + Ai2(m2 - cu^) 

+ (m - Up){mg + (M - LJp)/2){ujl + mg{M - Up)/2) 

+A?2(1 - m,(M - Up)/u^)/2] (3.6) 

as the exact counterpart of (2.43) under the same approximation. It is seen from (3.5) 
that the difference {n — p) is positive. 
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4 E.M. Mass Difference for the Nucleon 

The diagrams for the e.m. mass difference are given by figs. 3 (I, II, III) for a proton [uud) 
configuration to illustrate the underlying topology in accordance with the roles of the 
'active' and 'spectator' quarks in turn, as explained in Sec. 2. In each of these diagrams, 
two internal quark lines are joined by a photon line. The e.m. vertex at quark ^i has the 
strength e[l + 8x3 ]/6 from which the isospin matrix elements of a product of two such 
factors (shown for fig.3.111) have the forms 



< 



3rf)/6x(l + 3ri^V6|0';0"> (4.1) 



for the proton (uud) cnfiguration shown in III with ^3 as spectator, but in a basis (1;23) 
(which is consistent with the spin basis, eqs.(2.3-4)), corresponding to fig. la, viz. [37,39]: 

10' >= ■Ui(-U2'^3 - Usd2)/V2; 10" >= {-2diU2U3 + uid2U3 + uiU2d^)/\/Q (4.2) 

We note in parentheses that in fig. 3. Ill, the interchange of the two 'active' quarks 7^1 and 
7^2 does not give a new configuration, unlike in figs.l and 2; ((a) versus (b) configurations). 
It is now easy to check that the matrix elements <>' and <> " of (4.1) are 1/9 and 
— 1/9 for the proton configuration. After doing the corresponding neutron case, the two 
results may be combined in the single operator forms [39]: 

< . >'= (1 + 3r3)/36; <.>" = (!- 5r3)/36 (4.3) 

where r^ is the isospin operator for the nucleon as a whole [see eq.(3.5)], to be sandwiched 
between the neutron and proton states. The resultant isospin factor is then 

e^[< .>' + <.> "]/2 = 6^(1 - r3)/36 ^ -6^3/36 (4.4) 

After this book-keeping on the charge factors we can drop the isospin d.o.f. |0 > from 
the qqq wave function and, on the basis of the equality of the (.)' and (.)" contributions 
(2.22-23) for the spin matrix elements, it is enough to work with the (.)' type to represent 
the full effect. Collecting these details, the net isospin contribution to the e.m. (n — p) 
mass difference is just e^/18, which (of course) comes out with the correct (negative) sign 
in the resultant e.m. contribution to the total n — p difference after all the phase factors in 
the orbital-cum-spin space have been taken into account. The complete e.m. self energy 
of the nucleon (with operator t-^), with fig.3.III as the prototype, is now given by 

5M^ = 5:[-eV3/36]/(27r)^ f dWqud'qu'nnk-' x [23];[1];/(A3A2A^) (4.5) 

123 

where the various momentum symbols are as shown in fig. 3, with the primed quantities 
referring to the vertex on the right, but otherwise written in the same convention as in 
eqs.(2.5-10). The symbols within square brackets are analogous to (2.16-19): 

[ll = U{P')SAp[)t^,SF{p^)U{Py, {P' = P) (4.6) 

Tv 

[23];. = ^]^^[^"'75(M - i^.P'){mg - i^.p',)t^^{m, - t^.p2){M - ii.P)i,{m, + ii.p^)C] 

(4.7) 
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And the product of (4.6) and (4.7) works out as 

ME ^ ^!^|-^KAi + A;-A:^)(A2 + A',-P)/4 

-(Ai + a; - k^){mgUJ2 + rriqu'^ + 2cj2Cc'2)/2 

— (A2 + A2 — k'^){mqUJi + niqUj'i + 2ujiuj[)/2 

-(Ai + A2)(mg + uj[){mq + cj2)/2 - (A; + A'2)(mg + uoi){mq + cj2)/2 

-(A; + A2)(m, + cji)(m, + cj^)/2 - (Ai + A'2)(m, + uj[){mq + cj2)/2 

+ ("^g + lj{P -P^ - kY)[{mq + uJi){mq + uj^) + ("^g + cJi)(mg + 6^2) 

+ (mg + ooi){mq + 6^2) + {^q + cj^)(mg + cjg)]] (4.^ 



Some features of this "master" expression may be noted. There is a 'natural factor- 
ization'in the variables gi2 and q'12, except for the photon propagator /c~^, {k = qu — Qu)- 
Further, the two blobs are connected together by the 'spectator 'variable ^3 which is on 
the mass shell due to the presence of TgFg in eq.(4.3). 

The time-like (pole) integrations over each of d(Xi2 and 0-^2 can be carried out exactly a 
la (2.28) and its derivatives, since the 3D vertex function -Di20 in T^ does not involve au, 
etc. After this step qu, q'12 and ps are the 'right' 3D variables for the 'triple integration' 
whose essential logic may be stated as follows. The main strategy is to decouple the q 
and q' variables from the photon propagator k through the following device [28] : 

Since k is basically space- like, it is a good approximation to replace A;~^ by /c~^ which 
equals {qu — ^12)^? and drop the angular correlation in the two q- momenta (since the 
error in this neglect is zero in the first order [28]). Next we use the inequality [28] 

(a^ + b'^y^ < {2ab)-^; a -^ {qui, etc (4.9) 

which ensures the necessary factorizability in the q-variables.In principle the corrections 
to this inequality can be calculated since the neglected term is approximately equal to 
— (a — 6)^/(4a^6^) which is still factorizable, but this refinement is unnecessary in view of 
the smallness of the e.m. effect itself. After this simplification the rest of the integration 
procedure is straightforward since the q and q' integrations can be done analytically, 
and only a ID integration over \ps\ remains for numerical evaluation. The necessary 
expressions are collected in Appendix A and the numerical results for all contributions 
are given as under. 

The key parameters are the quark mass m^ and the size parameter (3"^, the latter being 
determined dynamically through the chain of eqs.(2.11-14). As noted in Seel already, 
the mass mq which is usually called the 'constituent' mass, should be viewed as the sum 
of the (flavour independent) 'mass function' m(p) for small p, plus a small "current mass" 
rric, in the spirit of Politzer additivity [22]. The mass function m(p) was generated in this 
BSE-cum-SDE framework through a Dynamical Chiral Symmetry Breaking mechanism 
in a non perturbative fashion [23]. Also from some related quark-loop calculations with 
qq mesons in recent times [24,28], it was found that for such 'low energy' processes the 
mass function m(p) is rather well approximated by m(0), so that [22], rriq = m(0) -|- rric. 
Therefore the d—u mass difference is the same at the 'constituent' or at the 'current' levels, 
and this is what has been denoted by 6m in the text (figs. 2). Its smallness compared to rUq 
justifies its neglect in all the functions except where it appears explicitly, viz., fig.2. We 
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take its value at 5m = 4MeV, as in related calculations [24,28], while the other quantities 
are predetermined from qq [14] and qqq [15] spectroscopy: 

mg = 2Q5MeV; f3\N) = 0.052GeV^ (4.10) 

so that there are no free parameters in the entire calculation. The results from Sees. 2-4 
are now summarized for (n-p) as : 

Ns^ = 5.5209 X IQ-^GeV-^^; [a = 1/3] (4.11) 

SMst = +l.7l3AMeV; 5M"< = -0A39QMeV. (4.12) 

Hence 

6M{net) = +1.28MeV; {vs.l.29MeV : Expt) (4.13) 

which is the principal result of this investigation, but subject to possible gauge invariance 
corrections [32]; see Appendix C and discussion in sec. 5. 3. 

5 Discussion, Summary and Conclusion 

This calculation fills up an important gap in the two-tier BSE formalism under 3D kernel 
support (CIA) for a simultaneous investigation of spectra and transition amplitudes of 
both qq and qqq varieties under a single umbrella [8,12]. 

5.1 Recapitulation 

To recapitulate the main points, the (first stage) 3D reductions of both the 2-body and 
3-body BSE's had yielded good agreement with the respective spectra [14,15], with a com- 
mon set of parameters Cq = 0.27 and ujq = 158Mel^ characterizing the non-perturbative 
gluon propagator, since the constituent mass ruq for spectroscopy [14,15], is essentially 
the dynamical mass function m{p) in the low momentum limit [22, 23]. 

More substantial tests of the formalism have come from the (second stage) reconstruc- 
tion of the 4D hadron-quark vertex function which carries the non-perturbative ofF-shell 
information in a closed form. This exercise was initially confined to the meson-gg vertex 
function whose exact reconstruction [8] had led to several useful results from 4D loop 
integrals for hadronic and e.m. transition amplitudes [8,16], to like integrals probing 
the momentum dependence of the quark mass function m{p) which is the 'chiral' limit 
(M^ = 0) [17,21,23] of the pion-quark vertex function. Indeed m{p) acts as the form 
factor for loop integrals determining the vacuum to vacuum transitions, and is found to 
predict correctly several condensates, from the basic < qq > [23] to 'induced' condensates 
[40] , under one roof. Further tests of the hadron-quark vertex function have come from 
SU(2) breaking effects like p — u mixing [24] and mass splittings in pseudoscalar mesons 
[28] , with only one additional parameter representing the d — u mass difference. 

The last link in our two-tier formalism has been a reconstruction of the 4D baryon-ggg 
vertex function on qq lines [8], to evaluate 3-body loop integrals. Although conjectured 
some time ago [12], a rigorous derivation [29] via Green's functions is outlined in Appendix 
B, in preparation for our main task: the n — p mass difference on identical lines to the 
qq case [28], viz., strong SU(2) breaking (fig 2) and e.m. contribution (fig 3), since the 
qqq vertex function is entirely determined by the same (gluon exchange) dynamics [7,23] 
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as qq. Strong SU(2) breaking is due to d — u mass difference, a la Weinberg [41], and 
more references on the physics of the problem may be found in [28]. The point to stress 
is that no free parameters are involved, so that the final value (4.13), although a single 
number, must not be treated as an isolated quantity, but as an integral part of a much 
bigger package. 

The 'QCD' status of this 3D-4D BSE formalism [8] viv-a-vis chiral perturbation theory 
[27] has already been explained in Seel: the gluon exchange character of the pairwise qq 
or qq interactions lends them a natural chiral invariance property at the input Lagrangian 
level with 'current' quarks. Thus the 'constituent' mass is not an input, but emerges as the 
low momentum limit of the dynamical mass function m{p) that characterizes the quark 
propagators appearing in the 2- and 3- body BSE's, as a result of DBxS [21, 17, 19, 23], 
since the 'current' masses of -u, d quarks give only a small additive contribution [22]. The 
empirical aspect of the gluon propagator concerns only its non-perturbative regime which 
often requires separate parametrization even in more orthodox BSE-SDE formulations 
[18]. In the present formulation, its explicit parametrization with two constants Cq and 
ujq) [23] symbolizes a 'closed form' representation of non-perturbative effects in the derived 
hadron-quark vertex function, but the returns are rich, especially the interlinkage of 4D 
loop integrals for different transition amplitudes [12,23,24,28], with the 3D BSE structures 
relevant to spectroscopy [14,15]. 

In contrast, Chiral Perturbation Theory [27] has an explicit QCD content, but re- 
lies more heavily on a perturbative treatment, as revealed by expansions in powers of 
small momenta and "current" masses rric [27] for a systemic derivation of the low energy 
structure of the Green's function in QCD [27]. It is a powerful method, highly success- 
ful in predicting items like ground state masses as well as their splittings, but its lack 
of a closed form representation prevents an equally successful prediction of 'soft' QCD 
effects in enough details, such as the momentum dependence of the mass function, or of 
hadron-quark vertex functions in general, and other observable effects such as L- excited 
spectra. 

5.2 Comparison with QCD-Sum Rules 

For a comparison with QCD-Sum Rules, it is useful to start by mentioning our neglect of 
the condensate contribution inserted in the internal quark lines, vide Fig Ic of [28], on the 
ground that it was found to be small in the 2-body case within the same BSE framework 
[28]. This contrasts sharply with the corresponding QCD-SR scenario [26] wherein the 
condensate contributions are the primary source of non-perturbative effects, as confirmed 
by explicit calculations [42]. Even more surprising is the inversion of the effect of the d — u 
mass difference on the n — p mass difference vis-a-vis the traditional low energy wisdom 
which requires rf — m > to make n — p > 0, as was the original motivation behind the 
famous Weinberg proposal [41]. (Incidentally our BSE result conforms to the Weinberg 
picture [41] and hence opposite of QCD-SR [42]). Indeed QCD-SR must rely heavily on 
the condensate contributions to compensate for the (negative) effect of the d — u mass 
difference, and leave a balance [42]. The e.m. effect too makes a comparable contribution 
to QCD-SR [42]. On the other hand the BSE approach, eq.(4.12), seems to predict only 
a modest e.m. effect which is about a fourth of, and of opposite sign to, the strong SU(2) 
breaking effect which dominates the entire scenario. 

What could account for such a sharp division between the two approaches? Perhaps 
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the reason should be sought in the very difference in the philosophy behind their respective 
premises: QCD-SR is basically a perturbative QCD approach designed from the high 
energy end, with the 'twist' terms (condensates) representing non-perturbative corrections 
of successively higher orders. In contrast the BSE is an intrinsically non-perturbative 
approach built from the low energy (spectroscopy) end, with the hard gluon exchange 
added perturbatively. The role of condensates in BSE is effectively subsumed in the 
constituent (non-perturbative) quark mass as well as the hadron-quark vertex function, 
so that any further condensate effects (fig Ic of [28]) in such a non-perturbative scenario 
can at best be residual [28]. This is not the case in QCD-SR wherein the quark condensates 
with their isospin splittings are the dominant source of non-perturbative effects on the 
n — p splitting [42]. The two scenarios are thus largely complementary, QCD-SR being 
rooted in hard QCD, and BSE-SDE in soft QCD premises respectively. 

5.3 Problem of Gauge Invariance of E.M. Effects 

Finally we come to a more vulnerable aspect of this investigation in company with the 
earlier study [28], viz., the lack of gauge invariance of the e.m. contribution (fig. (3)). 
Mercifully the e.m. contributions in both the meson [28] and baryon (present) cases are 
about a fourth of the u — d effect so as hopefully not to upset the overall stability of our 
result (4.13), but the need for a proper assessment of the g.i. corrections can hardly be 
overestimated. A general method for g.i. two-point functions for Qq systems has been 
given in [32], but it is not directly adaptable to extended vertex functions with arbitrary 
form factors, such as in the present situation. To handle these structures requires a 
different kind of strategy, a full-fledged formulation of which, is beyond the scope of 
this (already long) paper, and is best left to a separate communication. Nevertheless, 
as noted in Seel, we have made a beginning in this paper by outlining the main steps 
of the derivation of g.i. corrections for a typical two-body (kaon) case, in a simple and 
straightforward fashion, which amounts to the replacement of various momenta pi involved 
in the hadron-quark vertex functions by {pi — CiA), and expanding in powers of the e.m. 
fleld, with a view to calculate the additional diagrams a la [32]. Hopefully this method 
of generating e.m. gauge corrections is general enough to apply to other situations in 
which arbitrary momentum dependent hadron-quark vertex functions are involved. For 
the present situation of QED gauge corrections, the main steps are sketched in Appendix 
C, and the resultant correction to flg lb of [28] (i.e. flg la of [32]) is estimated to be 
—0.612MeV, which is 3/5 times, and of the same sign as, the e.m. value —1.032MeV 
for the mass difference K~ — K^ [28]. If this result is taken as a rough indication of the 
g.i. effect expected for the nucleon case on hand, it would mean a downward revision of 
the value (4.13) to about IMeV. However this is only a provisional estimate, pending a 
regular qqq calculation in the future. 

5.4 Conclusion 

To conclude, the principal motivation for this investigation, has been to demonstrate 
the practical feasibility of such realistic quark-loop calculations for the relativistic 3- 
quark problem within a full-fledged (BS) dynamical framework whose basic parameters 
are linked all the way to spectroscopy. The present calculation indeed suggests that not 
only quark-loops involving mesons [8,23,24,28] but even those involving the (less trivial) 

15 



qqq baryon are amenable to a similar degree of dynamical sophistication without excessive 
efforts, so that it makes sense to speak of an effective "4-fermion coupling" for both qq and 
qq pairs within a common parametric framework. This is somewhat remiscient of Bethe's 
"second principle" theory, originally suggested at the two-nucleon level of nuclear forces, 
now reinterpreted at the quark level, with a simple extension to include the antiquark 
in the dynamical description. (This extension would not make sense at the NN level, 
since the A^A^ and NN forces are very different from each other). Indeed such a dynamics 
had been strongly suggested (with concrete examples) in a perspective review not too 
long ago [43], but it seemed to have gone largely by default, as evidenced by a strong 
tendency in the contemporary literature to continue to rely on "ad-hoc form factors" [44] 
to simulate the vertex functions, instead of generating them dynamically. Hopefully, some 
efforts in this direction have been recently in evidence [45], using the Nambu Jonalasino 
model [21] of contact 4-fermion interactions (although such contact interaction models 
are probably too simple to realistically simulate confinement [29]). It is to be hoped that 
Bethe's "second principle" perspective will be upheld by such investigations, until such 
time as a fully satisfactory solution to QCD is forthcoming. 

The initial draft of this paper was prepared at the National Institute of Advanced 
Studies. We are grateful to Dr. Raja Ramanna for the warm NIAS hospitality. We also 
acknowledge Ms Chandana's help with some 'difficult' figures. One of us (ANM) is grate- 
ful to Prof. S.R.Chaudhury for a critical discussion on the g.i. problem in the e.m. 
contribution to the n — p mass difference. 

Appendix A: Evaluation of the Integral (4.5) 

The master expression (4.8) after being substituted in the full e.m. self energy contribution 
(4.5) is integrated over each dai2 and dai2'. The final result is 

^ 9 ^^i 2cu, 47r 47r An 2q,22qU 



123 



Fiqi2, qi2,p3)expi-'^[q!2 + & + ^PI]I P^) (A.l) 

where 

F{qi2,q'i2,p3) = 

{mg + uji2){mg + uju')k^ + (Mcug - -(M^ - ml))[{mg + ujuf 

+ {mq + uu'f + {nig + uji2){mg + Uu') - {nig + uJi2fDi2'/2uJi2' 
-{nig + u;i2')^^i2/2c^i2 - {rrig + uJi2){nig + Ui2')[Di2/2uJi2 
+Du'/2iuu']k^[{nig + 2LUi2){nig + 2iUu') - nil]/2 
-[{mg + 2uJi2){nig + 2uJu') - mJ][Di2/2cui2 + D^' /2loi2']/2 

k^-e[Du/uJi2 + Du'/LOi2] (A.2) 



lDuDi2' , ,4 ,2, 



8 uju^u' 

Using eq.(4.9), the integration over gi2 and g^2 can be done independently of each other, 
and thus can be written in a compact notation as follows 

2 



123 



U+^e^) I pI^P^^^^^^^F^ ^'-'^''^' (A-3) 
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where 

Fi = JnJu + [Mcu^ - -M^ + -ml]{2J2oJoo + 2JwJio) - 2J20I00 
— 2J10/10 + Ju Ju /2 — JoiJoiTng/2 — Jio Jio 
+mpooJoo - /00/00/2 - Jo2^o2/4 + /oi/oi (A.4) 

and with (n = 0,l,2;m = 0, 1, 2), 

J„^; Inm = 2-1/2 I ^^^^^^^ e(-i^-/'3')[V2gi2]"K + ^i2)"[l; ^^12/^12]; (A.5) 

J„J; /„^' = 2-^1^ j q,2dqi2 e(-i«-/^')[v^gi2]'"K + 2cc^i2)"[l; ^/^i2M2]; (A.6) 

Appendix B: Derivation of qqq Vertex Fn, Eq.(2.5) 

B.l: Method of Green's Functions 

We outhne here some essential steps leading to a formal derivation of eq.(2.5) which was 
written down in a semi-intuitive fashion in [12]. To that end we shall employ the method 
of Green's functions for 2- and 3- particle scattering near the hound state pole, since the 
inhomogeneous terms are not relevant for our purposes. For simplicity we shall consider 
identical spinless bosons, with pairwise BS kernels under CIA conditions [8] , first for the 
2-body case for calibration, and then for the 3-body system. 

B.2: Two-Quark Green's Function 

Apart from some results already goven in the text, we shall use the notation and phase 
conventions of [8,12] for the various quantities (momenta, propagators, etc). The 4D 
qq Green's function G{piP2;piP2') near a bound state satisfies a 4D BSE without the 
inhomogeneous term, viz. [8,12], 

i(2nYG{piP2;pi'p^) = ^i-'/^2-' j dp^'dp^'K{p^P2;pi"p2")G{pi"p^'-pi'p2') (B.2.1) 

where 

Ai=Pi' + V, (B.2.2) 

and iTLq is the mass of each quark. Now using the relative 4- momentum q = {pi — P2)/2 
and total 4-momentum P = pi + P2 (similarly for the other sets), and removing a 5- 
function for overall 4-momentum conservation, from each of the G- and K- functions, 
eq.(B.2.1) reduces to the simpler form 

i{2TxfG{q.q') = Ar^ A2-' / dq"Mda"K{q, q")G{q\ q') (B.2.3) 

where q^ = q^ — crPfj,, with a = [q.P)/P'^, is effectively 3D in content (being orthog- 
onal to P^). Here we have incorporated the ansatz of a 3D support for the kernel K 
(independent of a and a'), and broken up the 4D measure dq" arising from (2.1) into the 
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product dq"Mda" of a 3D and a ID measure respectively. We have also suppressed the 
4-monientum P^ label, with (P^ = — M^), in the notation for G{q.q'). 
Now define the fully 3D Green's function G{q, q') as 



G{q,q') = M^dada'G{q,q') (B.2.4) 
and two (hybrid) 3D-4D Green's functions G{q,q'), G{q,q') as 

G{q, q') = J MdaG{q, q'); G{q, q') = J Mda'G{q, q'); (B.2.5) 
Next, use (B.2.5) in (B.2.3) to give 

t{27r)'G{q, q') = Ar^Aa"^ / dq"K{q, q")G{q" , q') (B.2.6) 
Now integrate both sides of (B.2.3) w.r.t. Mda and use the result [8] 

[ MdaAi-^A2'^ = 27riD-\qy, D{q) = 4uj{u;^ - M^/A); uj^ = rrig'^ + f (B.2.7) 
to give a 3D BSE w.r.t. the variable q, while keeping the other variable q' in a 4D form: 

{2nf G{q, q') = D'' J dq"K{q, q")G{q", q') (B.2.8) 

Now a comparison of (B.2.3) with (B.2.8) gives the desired connection between the full 
4D G-function and the hybrid G(g, g')-function: 

2mG{q, q') = D(g) Ar' A2-^G'(g, q') (B.2.9) 

Again, the symmetry of the left hand side of (B.2.9) w.r.t. q and q' allows us to write the 
right hand side with the roles of q and q' interchanged. This gives the dual form 

2mG{q, q') = D{q')/\r' A^'Ciq, q') (B.2.10) 

which on integrating both sides w.r.t. Mda gives 

2mG{q, q') = D{q')Ar'Ar'G{q, q'). (B.2.11) 

Substitution of (B.2.11) in (B.2.9) then gives the symmetrical form 

{2mfG{q, q') = D{q)A^-'A2-'G{q, q')D{q')Ar'Ar' (B.2.12) 

Finally, integrating both sides of (B.2.8) w.r.t. Mda', we obtain a fully reduced 3D BSE 
for the 3D Green's function: 

i2nfGiq,q') = D'^q j dq" K{q,q")G{q\q') (B.2.13) 

Eq. (B.2.12) which is valid near the bound state pole (since the inhomogeneous term has 
been dropped for simplicity) expresses the desired connection between the 3D and 4D 
forms of the Green's functions; and eq(B.2.13) is the determining equation for the 3D 



form. A spectral analysis can now be made for either of the 3D or 4D Green's functions 
in the standard manner, viz., 

G(g, q') = y: Mq; P)Kiq'; p)/{p' + m') (b.2.i4) 

n 

where $ is the 4D BS wave function. A similar expansion holds for the 3D G-function 
G in terms of (pniQ)- Substituting these expansions in (B.2.12), one immediately sees the 
connection between the 3D and 4D wave functions in the form: 

2m^q, P) = Ar'A2-^/^(g)0(g) (B.2.15) 

whence the BS vertex function becomes T = D x (f)/{2ni) as found in [8]. We shall make 
free use of these results, taken as qq subsystems, for our study of the qqq G-functions in 
Sections 3 and 4. 

B.3: 3D Reduction of the BSE for 3-Quark G-function 

As in the two-body case, and in an obvious notation for various 4-momenta (without 
the Greek suffixes), we consider the most general Green's function G{pip2P3;pi'p2'p3') for 
3-quark scattering near the bound state pole (for simplicity) which allows us to drop the 
various inhomogeneous terms from the beginning. Again we take out an overall delta 
function S{pi + P2 + Ps — P) from the G-function and work with two internal 4-momenta 
for each of the initial and final states defined as follows [12]: 

^3^3 =pi-p2 ; Sr]3 = -2p3 +P1+P2 (B.3.1) 

P = Pl+P2+P3=Pl +P2 +P3 (B.3. 2) 

and two other sets C,i,tii and ^2, V2 defined by cyclic permutations from (B.3.1). Further, as 
we shall consider pairwise kernels with 3D support, we define the effectively 3D momenta 
Pi, as well as the three (cychc) sets of internal momenta ^i,r/i, (i = 1,2,3) by [12]: 

Pi=Pi- i^iP ; ii = ^i- SiP ; fji - tiP (B.3. 3) 

z/, = {P.pi)/P^ ; s, = {P.^i)/P' ; t, = {P.v^)/P' (B.3.4) 

v3s3 = ui — 1^2 ', 3t3 = — 2z/3 + z/i -|- z/2 (-|-cyclicpermutations) (B.3. 5) 

The space-like momenta pi and the time-like ones z/j satisfy [12] 

Pi+P2+P3 = ; z/i + z/2 + z/3 = 1 (B.3. 6) 

Strictly speaking, in the spirit of covariant instantaneity, we should have taken the relative 
3D momenta ^,fj to be in the instantaneous frames of the concerned pairs, i.e., w.r.t. the 
rest frames of Pij = Pi + Pj ; however the difference between the rest frames of P and Pij 
is small and calculable [12], while the use of a common 3-body rest frame (P = 0) lends 
considerable simplicity and elegance to the formalism. 

We may now use the foregoing considerations to write down the BSE for the 6-point 
Green's function in terms of relative momenta, on closely parallel lines to the 2-body case. 
To that end note that the 2-body relative momenta are qij = {pi —pj)/2 = \^C,k/2, where 
(ijk) are cyclic permutations of (123). Then for the reduced qqq Green's function, when 
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the last interaction was in the (ij) pair, we may use the notation G{^kVk] C,kVk'), together 
with 'hat' notations on these 4- momenta when the corresponding time-hke components 
are integrated out. Further, since the pair C,k,Vk is permutation invariant as a whole, 
we may choose to drop the index notation from the complete G-function to emphasize 
this symmetry as and when needed. The G-function for the qqq system satisfies, in the 
neighbourhood of the bound state pole, the following (homogeneous) 4D BSE for pairwise 
qq kernels with 3D support: 

ti27r)'Gi^r]; ^V) = E Ar^A^"^ / dq'l,Mdau"K{q,2, q'l2)G{C3m"; ^3^') (B.3.7) 

123 

where we have employed a mixed notation (gi2 versus ^3) to stress the two-body nature 
of the interaction with one spectator at a time, in a normalization directly comparable 
with eq.(B.2.3) for the corresponding two-body problem. Note also the connections 



0'12 



V3S3/2; q,2 = Vsis/^; m = -P3, etc (B.3.8) 



The next task is to reduce the 4D BSE (B.3.7) to a fully 3D form through a sequence 
of integrations w.r.t. the time-like momenta Sj,tj applied to the different terms on the 
right hand side, provided both variables are simultaneously permuted. We now define the 
following fully 3D as well as mixed (hybrid) 3D-4D G-functions according as one or more 
of the time-like C,,f] variables are integrated out: 



Gi^fj; i'f,') =J J J J dsdtds'dt'Gi^f]- eV) (B.3.9) 

which is S'3-symmetric. 

Gsni^mm =J Jdtsdts'Gi^r^-^W); (B.3.10) 

G3div;iW) = I Jdssdss'Gi^r^-er]'); (B.3.11) 

The last two equations are however not symmetric w.r.t. the permutation group S3, since 
both the variables C,, tj are not simultaneously transformed; this fact has been indicated in 
eqs.(B.3.10-ll) by the suffix "3" on the corresponding (hybrid) G-functions, to emphasize 
that the 'asymmetry' is w.r.t. the index "3". We shall term such quantities "5'3-indexed" , 
to distinguish them from 5'3-symmetric quantities as in eq. (B.3.9). The full 3D BSE for the 
G-function is obtained by integrating out both sides of (B.3.7) w.r.t. the st-pair variables 
dsidsj'dtidtj' (giving rise to an S'3-symmetric quantity), and using (B.3.9) together with 
(B.3.8) as follows: 

(27r)=^G(er); ^'r)') = ^ ^"'(^12) / dq';2K{q^2, q'U)G{i"f^"- i'f^') (B.3.12) 

123 

This integral equation for G which is the 3-body counterpart of (B.2.13) for a qq system 
in the neighbourhood of the bound state pole, is the desired 3D BSE for the qqq system in 
a. fully connected iorm, i.e., free from delta functions. Now using a spectral decomposition 
forG 

G{iv; i'v) = E Mim P)m'v; p)/{p' + m') (B.3.13) 
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on both sides of (B.3.12) and equating the residues near a given pole P^ = — M^, gives 
the desired equation for the 3D wave function for the bound state in the connected 
form: 

(271)' m;P) = J2D-\qi2) J dq';,Kiqu,qum"v";P) (B.3.14) 



123 



Now the S's-symmetry of in the {C,i,'i]i) pair is a very useful result for both the solution 
of (B.3.14) and for the reconstruction of the 4D BS wave function in terms of the 3D wave 
function (B.3.14), as is done in the subsection below. 

B.4: Reconstruction of the 4D BS Wave Function 

We now attempt to re-express the 4D G-function given by (B.3.7) in terms of the 3D G- 
function given by (B.3.12), as the qqq counterpart of the qq results (B. 2. 12-13). To that 
end we adapt the result (B.2.12) to the hybrid Green's function of the (12) subsystem given 
by Gsr], eq.(B.3.10), in which the 3-momenta fj^jfj'^ play a parametric role reflecting the 
spectator status of quark ^^3, while the active roles are played by qu, qu' = V^i^s, ■CsO/^, 
for which the analysis of subsec.B.2 applies directly. This gives 

{2m)'Gs,{^sV3;^3%) = D{q^2)Ar'A2-'G{is7^s;i3%)D{q[,)Ar'Ar' (B.4.1) 

where on the right hand side, the 'hatted' G-function has full S'3-symmetry, although (for 
purposes of book-keeping) we have not shown this fact explicitly by deleting the suffix 
'3' from its arguments. A second relation of this kind may be obtained from (B.3.7) by 
noting that the 3 terms on its right hand side may be expressed in terms of the hybrid 
G^ functions vide their definitions (B.3.11), together with the 2-body interconnection 
between (^3,^3') and (^3,^3) expressed once again via (B.4.1), but without the 'hats' on 
?73 and 773'. This gives 

{V3mfG{CsV3;^3m) = iVSmyGi^r];ev') 

= J2^^''^2-\mV3) fdq';,Mdau"K{qu,q';2)Gi^s"v/;aV3) 

123 "^ 

= Y. D{q^2)AL-'A2-'Gsdi3V3; i'sV3')Ar'Ar' (B.4.2) 

123 

where the second form exploits the symmetry between C,,ti and ^',77'. 

At this stage, unlike the 2-body case, the reconstruction of the 4D Green's function 
is not yet complete for the 3-body case, as eq. (B.4.2) clearly shows. This is due to the 
truncation of Hilbert space implied in the ansatz of 3D support to the pairwise BSE kernel 
K which, while facilitating a 4D to 3D BSE reduction without extra charge, does not have 
the complete information to permit the reverse transition (3D to 4D) without additional 
assumptions; see [29] for details. The physical reasons for the 3D ansatz for the BSE 
kernel have been discussed in detail elsewhere [23,29], vis-a-vis contemporary approaches. 
Here we look upon this "inverse" problem as a purely mathematical one. 

We must now look for a suitable ansatz for the quantity G^^ on the right hand side 
of (B.4.2) in terms of known quantities, so that the reconstructed 4D G-function satisfies 
the 3D equation (B.3.12) exactly, as a "check-point" for the entire exercise. We therefore 
seek a structure of the form 

G^ii3V3; iW) = Giism; i'^v's) x ^(^3,^3') (B.4.3) 
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where the unknown function F must involve only the momentum of the spectator quark 
7^3. A part of the tj^iTj^' dependence has been absorbed in the G function on the right, 
so as to satisfy the requirements of S'3-symmetry for this 3D quantity [29] . 

As to the remaining factor F, it is necessary to choose its form in a careful manner so 
as to conform to the conservation of 4-momentum for the free propagation of the spectator 
between two neighbouring vertices, consistently with the symmetry between p^ and p^' . 
A possible choice consistent with these conditions is the form (see [29] for details): 

^(P3, Ps') = C3 A3-^5(z/3 - y^') (B.4.4) 

Here A3~^ represents the "free" propagation of quark ^^3 between successive vertices, 
while C3 represents some residual effects which may at most depend on the 3-momentum 
P3, but must satisfy the main constraint that the 3D BSE, (B.3.12), be exp/«cz% satisfied. 
To check the self-consistency of the ansatz (B.4.4), integrate both sides of (B.4.2) w.r.t. 
ds^ds^' dt^dt^' to recover the 3D S'3-invariant G-function on the left hand side. Next, in 
the first form on the right hand side, integrate w.r.t. ds^ds^' on the G-function which 
alone involves these variables. This yields the quantity G^^. At this stage, employ the 
ansatz (B.4.4) to integrate over dt^dt^'. Consistency with the 3D BSE, eq. (B.3.12), now 
demands 

C3 I I du^du^^^-^5{u^ - z/3') = 1; {sincedt = du) (B.4.5) 



The ID integration w.r.t. dv^ may be evaluated as a contour integral over the propagator 
A"-*^ , which gives the pole at z/3 = 0)3 /M, (see below for its definition). Evaluating the 
residue then gives 

C3 = m/{MCj^)- CjI = niq^ + p2 (B.4.6) 

which will reproduce the 3D BSE, eq. (B.3.12), exactly] Substitution of (B.4.4) in the 
second form of (B.4.2) finally gives the desired 3-body generalization of (B.2.12) in the 
form 

3G(er7;eV) = Y.D{qi2)^iF^2FD{c^,^)^iF' ^2F'G{i^f,^-A'A)[^^Fl{Mnu^)] (B.4.7) 

123 

where for each index, A^r = — iA"^ is the Feynman propagator. 

To find the effect of the ansatz (B.4.4) on the 4D BS wave function ^{^1]; P), we do 
a spectral reduction like (B.3.13) for the 4D Green's function G on the left hand side of 
(B.4.2). Equating the residues on both sides gives the desired 4D-3D connection between 
$ and 0: 

H^V, P) = E /^(gi2) A,- ArV(e^; P) X \^^^^^^ (B.4.8) 

^From (B.4.8) and eq.(2.1) of the text, we infer the structure of the baryon-ggg vertex 
function V3 as given in eq.(2.5) of the text. For a detailed discussion of the significance 
of this result, vis-a-vis contemporary approaches, see [29]. 

Appendix C: Gauge Corrections to Kaon E.M. Mass 

We outline here a practical procedure to evaluate the gauge corrections to the e.m. self- 
energy of a qq system, vide fig. lb of [28], pending a more systematic treatment in a 
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later paper. This two-body exercise should hopefully serve as a fascimile of the effect 
expected for the present qqq case. For brevity we shall refer to the figures of KL [32] 
in their notation without drawing them anew. Thus fig lb of [28] corresponds to fig la 
of KL [32], except for the presence of the hadron lines at the two ends. We shall call 
this simply 'la', with the understanding that the hadron lines are 'attached' to la. For 
the actual mathematical symbols (including phase conventions) we shall draw freely from 
[28], without explanation. In [28], only la of [32] was calculated, but now one must add 
2(a,b,c,d,e) of [32], all with hadron lines understood at the two ends of each. There is 
no need to calculate lb or Ic of [32] which are mere e.m. self-energies of single quarks 
(g.i. by themselves), and are routinely absorbed in quark mass renormalization (of little 
significance in this phenomenological study which has these masses as inputs). 

A new ingredient is a 4-point vertex in each of 2(a,b,c,d), and two 4-point vertices in 
2e, except that the word 'point' is now understood as an extended structure characterized 
by the hadron-quark vertex function D{q)(j){q) where one must insert a photon line in 
each such Hqq blob. Since it is not a standard point vertex, the method [32] of inserting 
exponential phase integrals with each current is not technically feasible; instead we may 
resort to the simple-minded substition pi — eiA{xi) for each 4-momentum pi (in a mixed 
p, X representation) occurring in the structure of the vertex function, which has the same 
physical content, at least up to first order in the e.m. field, without further comment. 
This amounts to replacing each q^ occurring in T{q) = D{q)(f){q), by q^ — CgA^, where 
Cq = 771261 — 77^162. The net result in the first order in A^ is a first order correction to 
r(g) of amount eqj{q).A defined by 

j{q).A = -4Mg.A0(g)(l - D{q)/{4M(3^)) (C.l) 

where we have made free use of various symbols and definitions in [28]. (The effect of 
the hat structure of q on the e.m. substitution is ignored in this approximate treatment). 
This effective 4-point vertex function is operative at one end in each of 2a,2b,2c,2d of KL 
[32] and at both ends of 2e. For the e.m. vertex at the quark lines of 2(a,b,c,d), we use 
simply iCi'-f.A, as in [28]. The matrix elements can now be written down on exactly the 
same lines, and the same phase convention as in [28] to keep proper track of the gauge 
corrections with sign. We need write these down only for 2a and 2e, noting the equalities 
2a=2b, as also 2c=2d, and the further substitutions (1) -^ (2) and vice versa to generate 
2c(=2d) from 2a(=2b). The contribution from 2a [32] to the e.m. quadratic self-energy 
of a kaon is expressible as 

Ml = Nl{2Ti)-'e,eq j j{q)^D{q')(t>{q')Tr[^,DF^,{k) 

Sf{pi - rhik)iei-iySF{Pi)l^^SF{-P2)]d^qd'^k (C.2) 

where p'^ = pi+rh2k and P2 = p'2 = P2 — '^2k are the 4-momenta of the quarks at the other 
(right-hand) end, and the photon propagator in the Landau gauge is —i{6^u — kfj,ki,/k'^)/k'^. 
To make better use of the techniques outlined in [28], it is convenient to change the variable 
from k^ to q'^, noting that q' = q + m2k, which gives d'^k = d'^q'/rh^, etc. This shows 
that fig 2a(=2b), where the photon line ends on the heavier quark mi, gives a bigger 
contribution than does fig.2c(=2d) which would give rh{^ arising from the d'^k-measnie. 
Evaluating the traces, and integrating over the poles of the two time-like momenta go 
and go gives for the sum of the contributions from 2a-2d to the quadratic mass difference 
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between K^ and K as a product of two 3D quadratures after some simplifications with 
factorable approximations a la [28]: 

QNlM5{e,e,) r ^. t-., 00' ^^ D{q) ^ 



5M. 



2 






[(g'(2 - 4/7r) - qq'/3){M^ - 5m' + D{q')cu['' /2 + D{q')uj'^'^ /2) 
+^m2qq'iD{q')iu'f^/2 + M^ - 5m')] + [1^2] (C.3) 

Here 5{eieq) is the K^ minus K~ difference between the indicated charge factors associated 
with line 'i', while cufg = ""^i 2 + ^'^ ^.nd ujfj^ = m\ + {q — rhik)'. 

Next the contribution to 5M' arising from fig 2e of KL [32] which involves the product 
of two vertex blobs like (CI) is given by 

5Ml = iNl{2Tiy'>e,' J d''qd*kDF^,{k)j{q)^j{q),Tr[^,SF{p^ - m^k)-f,SF{-p2 + m^k)] 

(C.4) 
This integral is somewhat different in structure from (C.2) in as much as /c^ is fully 
decoupled from either wave function 0, 0', both of which have the same argument q. This 
makes it possible to integrate first over d'^k as well as the time-like component go of g^ 
neither of which is involved in the vertex function. The relevant integral after tracing and 
rearranging has the form 

F{q) = 3{-t)' Jd'kJdqok~'{5^,-k^k,/e) 

[f -ql + mim2 - mim2{P - kf]/{AiA2) (C.5) 

where Aj = mf + {pi — rhik)'. The integral which is entirely convergent works out after 
some standard manipulations involving Feynman techniques as well as differentiation 
under integral signs as 

F{q) = 67r^[mim2 + f + A] [VX - ^A - mim2M^]/{mim2M)' (C.6) 

where A = mim2M' + D{q)/2M. And the final expression for (C.4) in terms of (C.6) is 

5Ml = Nl{27^)-H{el) fd^m'Fiq) (C.7) 



Further evaluation of (C.3) and (C.7) can be made a la [28] in a straightforward way. 
The key ingredients are 

5eieg = 0.236e'; (5e2eg = 0.1396^; (JeJ = -0.0294el (C.8) 

The break-up of the final results for the diagrams 2(a-e) after dividing the results of(C.3) 
and (C.7) by 2M, since 5M' = 2M5M, is (in MeV): 

(^M2a+26 = -0.6996; (5M2c+2d = +0.1358; 5M2e = -0.0481; 5 M^ = -0.612 MeV. 

(C.9) 
All these corrections, which reinforce one another due to a complex interplay of signs, add 
up to a figure which increases the value -1.032 MeV due to fig 1(b) found in [28], to -1.644 
MeV. A more comprehensive paper for the gauge corrections to fig 1(a) of [32] for the 
other mesons, as well as the qqq baryon, will be communicated separately. This estimate 
has been used, on a percentage basis, as a facsimile for the gauge correction expected for 
the e.m. self energy of the nucleon, in sec. 5. 3 of text. 
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Figure Captions 

Fig.l: Diagrams for BS normalization of Baxyon-qqq vertex function. 1(a) shows quark 
#1 emitting a zero momentum photon (/c = 0); its last qq interaction was with ^2, while 
7^3 is the spectator. 1(b) is the same diagram with the roles of 4^1 and #2 interchanged. 
1(c) denotes schematically two more such pairs of diagrams obtained with cyclical per- 
mutations of the indices (123) in pairs. The 4- momenta on the quark lines are shown as 
used in the text. 

Fig. 2: Diagrams for the two-point interactions of the quark lines with the mass shift 
operator —5mT^^^y2 in place of the photon in fig.l, but otherwise with identical topolog- 
ical correspondence of figs.2(a,b,c) to figs.l(a,b,c). 

Fig. 3: Diagrams for the e.m. self-energy of the uud (proton) configuration. 3(111) is 
shown in detail with full momentum markings as employed in the text, and corresponds 
to quark #3 as the spectator, while the quark lines #1 and #2 are joined by a transverse 
photon line. Similarly 3(1) and 3(11) correspond to #1 and #2 respectively as spectators 
in turn. Note that, unlike in fig.l and fig. 2, the interchange of #1 and #2 in fig.3(III) 
does not give a new configuration. 

REFERENCES 

[1] L.D.Faddeev, Sov.Phys.JETP 12, 1014 (1961). 

[2] D. Z. Freedman, et al, Nuovo Cim.43A, 258 (1966). 

[3] R.Blankenbecler, et al, Nucl.Phys.l2, 629 (1959). 

[4] R.P.Feynman et al, Phys.Rev.D3, 2706 (1971). 

[5] D.B.Lichtenberg, Phys.Rev.178, 2197 (1969). 

[6] a)M.Levy, Phys.Rev.88,72 (1952); b)A.Logunov and A.N.Tavkhelidze, Nuovo Ci- 
mento 29, 380 (1963); c)R.Blankenbecler and R.Sugar, Phys.Rev.142, 105 (1966); 
d)V.Kadychevsky, Nucl.Phys.B6, 125 (1968). 

[7] A.N.Mitra, Zeits.f.Phys.C8, 25 (1981). 

[8] A.N.Mitra and S.Bhatnagar, Int.J.Mod.Phys.AT, 121 (1992). 

[9] Yu.L.Kalinowski et al, Phys.Lett.B231, 288 (1989). 

[10] M.A.markov, Sov.J.Phys.3, 452 (1940); 
H.Yukawa, Phys.Rev.77, 219 (1950). 

[11] Particle Data Group, Phys.Rev.D54 July 1-Part I (1994). 

[12] A.N.Mitra and I.Santhanam, Few-Body Syst.l2, 41 (1992). 

[13] R.Barbieri and E. Rimmidi, Nucl.Phys.B141, 413 (1978). 

[14] K.K. Gupta et al, Phys.Rev.D42, 1604 (1990). 



25 



[15] A. Sharma et al, Phys.Rev.D50, 454 (1994). 

[16] N.N.Singh et al, Phys.Rev.D38, 1454 (1988); 

A.N.Mitra, Proc.Ind.Nat.Sci.Acad.54A, 179 (1988). 

[17] S.Adler and A.C.Davies, Nucl.Phys.B244, 469 (1984). 

[18] For a review of QCD motivated SDE-BSE, see C.D.Roberts and A. G. Williams, 
Prog.Part.Nucl.Phys.31, 477-575 (1994). 

[19] A.Le Youanc et al, Phys.Rev.D29, 1233 (1984). 

[20] A.Mittal and A.N.Mitra, Phys.Rev.Lett.57, 290 (1986). 

[21] Y.Nambu and G.Jona-Lasino, Phys.Rev.121, 345 (1961). 

[22] 22) H.D. Politzer, Nucl.Phys.B117, 397 (1976). 

[23] A.N.Mitra and B.M.Sodermark, Int.J.Mod.Phys.A9, 915 (1994). 

[24] A.N.Mitra and K.-C. Yang, Phys.Rev.C51, 3404 (1995). 

[25] W.BuechmuUer, Ed. Current Physics Sources and Comments 9 Quarkonia, North- 
Holland 1992. 

[26] V.L.Chernyak and A.R.Zitnitsky, Nucl.Phys.B246, 52 (1984); 
B.L.Ioffe and A.V.Smigla, Nucl.Phys.B232, 109 (1984). 

[27] I.Gasser and H.Leutwyler, Ann.Phys.(N.Y.)158, 142 (1984); 

J.Gasser and H.Leutwyler, Nucl.Phys.B250, 465, 517, 539 (1985). 

[28] A.N.Mitra, Intl.J.Mod.Phys.All, 5245 (1996). 



[29] A.N.Mitra, LANL |hep-th/9803U^ : 3D-4D Interlinkage oi qqq Wave Functions Under 



3D BS Kernel Support. 

[30] S.Weinberg, Phys.Rev.133B, 232 (1964). 

[31] J.M.Blatt and V.F.Weisskopf, Theoretical Nuclear Physics, John Wiley and Sons Inc, 
New York; Chap 2(C) (1952). 

[32] L.Kisslinger and Z.Li, Phys.Rev.Lett.74, 2188 (1995). 

[33] A.N. Mitra and A.Mittal, Phys.Rev.D29, 1399 (1984). 

[34] A.N.Mitra, A.Sharma and B.M. Sodermark, Few-Body Syst.l9, 1(1995); 

[35] J.Bijtebier, Nuovo Cimento 81 A, 423 (1985). 

[36] M. Verde, Handbuch der Physik, 39, 170 (1957). 

[37] A.N. Mitra and M.H. Ross, Phys.Rev.158, 1630 (1967). 

[38] K.Nishijima, A.S. Singh, Phys.Rev.1621, 1740 (1967); 
I.Santhanam et al, Int.J.Mod.Phys.E2, 219 (1993). 

26 



[39] A.N. Mitra and S. Sood, Fortschritte der Physik 25, 649-744 (1977). 

[40] W.Y.P.Hwang et al, Chinese J. of Phys.33, 397 (1995). 

[41] S.Weinberg, I.Rabi Festschrift, ed. L.Motz (New York Academy of Sciences), pl85 
(1977). 

[42] K.-C.Yang, Topics in QCD sum Rules, Ph.D. Thesis, Nat Taiwan Univ 1994. 

[43] I.Santhanam et al, Few-Body Syst.7, 1 (1990). 

[44] M.A.Ivanov et al, Few-Body Syst.21, 131 (1996); also for several other references. 

[45] S.Huang and J.Tjon, Phys.Rev.C49, 1702 (1994); 
N.Ishii et al, Austr J Phys 50, 123 (1997); 
W.Bentz, J. Korean Phys. Soc.29 Suppl, 5352 (1996). 



27 



